Entanglement dynamics in the Lipkin-Meshkov-Glick model 
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The dynamics of the one-tangle and the concurrence is analyzed in the Lipkin-Meshkov-Glick 
model which describes many physical systems such as the two-mode Bose-Einstein condensates. We 
consider two different initial states which are physically relevant and show that their entanglement 
dynamics are very different. A semiclassical analysis is used to compute the one-tangle which 
measures the entanglement of one spin with all the others, whereas the frozen-spin approximation 
allows us to compute the concurrence using its mapping onto the spin squeezing parameter. 
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I. INTRODUCTION 

In the recent years, the interplay between entangle- 
ment and quantum phase transitions has been the subject 
of many studies. Since the pioneeringworks on the Ising 
model in a transverse magnetic field 1 , 2] exhibiting the 
key role played by entanglement in quantum critical phe- 
nomena, lots of effort have been devoted to characterize 
the intricate structure of the ground state in spin systems 
J441A4 iL 1, 1, 10, li, 12, 13, ii, 15, 16.^ 17^8, li 
12(1 I21L |2^ as well as in electron models |2J, |2J, [23 ■ AH 
these studies show that entanglement, a genuine quan- 
tum property, can be used to probe the phase diagram of 
a system and especially to detect quantum phase transi- 
tions. Unfortunately, such analyses are often restricted 
to one-dimensional (ID) systems where exact solutions 
exist and allow one to deal with the thermodynamical 
limit. Usually, higher-dimensional systems require a nu- 
merical treatment which is limited by the Hilbert space 
dimension which grows exponentially with the number 
of degrees of freedom. However, when the Hamiltonian 
has a large symmetry group as is here the case, it is still 
possible to study the large- iV limit. 

In the present paper, we consider N mutually interact- 
ing spins 1/2 embedded in a magnetic field. The restric- 
tion of this model to the fully symmetric subspace can be 
mapped, via the Schwinger representation, onto interact- 
ing bosons in a two-level system coupled by a tunneling 
term. Thus, although it has initially been introduced 
by Lipkin, Meshkov, and Ghck (LMG) 26, 27, 28] in 
nuclear physics, this model is also relevant to describe 
the Josephson effect, or two-mode Bose-Einstein conden- 
sate (BEC). This ubiquity in very different domains is 
certainly one of the reason why it has been periodically 
rediscovered (see, for example, [23, Hfl IHl ) . Its integra- 
bility has even recently been proved in a series of paper 
(see Ref. 62\ for a review) and a complete solution has 
been derived using the algebraic Bethe ansatz [33 ■ The 
ground state entanglement properties of this model have 
been analyzed for both ferromagnetic |34L l35| and an- 
tiferromagnetic coupling |3y| using the concurrence |33 
(see below for details). Here, our goal is to analyze the 
entanglement dynamics which has, so far, been mostly 



investigated in ID spin systems |38ll39|. 

The possibility of experiments in BEC's has led sev- 
eral groups to study the dynamics in the LMG model 
but ent ang lement has been discussed only recently |43, 
I41I 143 . [43 . In this study, we investigate the time evo- 
lution of two peculiar states which are relevant in the 
BEC context since they correspond to situations where 
all atoms are either in one of the two modes or equally 
distributed between both modes. For these two states, 
we discuss the entanglement dynamics through two ob- 
servables: namely, the one-tangle which measures the en- 
tanglement of one spin with all others and the concur- 
rence which quantifies the two-spin entanglement. In the 
next section, we introduce the LMG model and discuss 
its phase diagram. In Sec. IIIIL we define the entangle- 
ment measures used throughout this paper and give sev- 
eral canonical examples. The ground-state entanglement 
properties |3J,|33,|3a| are briefly recalled in Sec. IIVI The 
quantum dynamics is presented in Sec. for the two 
states mentioned above. We show that the one-tangle 
and the concurrence have completely different behaviors 
depending on the initial state. These characteristics can 
be understood in several limiting cases that we discuss in 
detail. 



II. THE LIPKIN-MESHKOV-GLICK MODEL 

We consider a system of mutually interacting spins 1/2 
embedded in a magnetic field described by the following 
Hamiltonian introduced by LMG [23, 1271 128| in a more 
general form: 



H 



2A 



hJ2<'l 



— ~M^ ^hSz + — , 



(1) 
(2) 



where the CTq's are the Pauli matrices and Sa — X^i ^a/"^- 
The prefactor 1/N is necessary to get a finite free energy 
per spin in the thermodynamical limit. The Hamiltonian 
H preserves the magnitude of the total spin and does not 
couple states having a different parity of the number of 



spins pointing in the magnetic field direction (spin-flip 
symmetry): namely, 



[H,S'] = 0, 



^H' 



(3) 
(4) 



The phase diagram of this model can be easily derived 
using a semiclassical description |30j . For a ferromagnetic 
couphng (A > 0), the system undergoes a second-order 
quantum phase transition at the critical field he — A, 
whereas for antiferromagnetic interactions, a first-order 
quantum phase transition occurs at zero field. 



III. ENTANGLEMENT MEASURES 



Note that for pure two-spin states, we have an exact 
equality. In the present study, we have systematically 
checked that this conjecture ||SJ) was always verified. 

Finally, we would like to stress that tracing out the 
density matrix p can induce or destroy some correlations 
between spins. For instance, a nonseparable state, such 
as the Greenberger-Horne-Zeilinger (GHZ) state jig 

|GHZ) = |^) = -^(|TTT> + UU)), (9) 

has a maximum one-tangle (t — 1) and a vanishing con- 
currence (C = 0). It is thus important to keep in mind 
that these measures only give partial information about 
the entanglement of the state. 



To analyze the entanglement dynamics, we focus here 
on two commonly used quantities: the one-tangle which 
measures the entanglement of one spin with the others 
and the concurrence which measures the two-spin entan- 
glement. For a given pure state {ip), both quantities are 
computed from the density matrix p = \'ip){ip\. 

The one-tangle is defined as follows: 



4 det p. 



(1) 



1 
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.(1) 



where pl is the one-spin reduced density matrix ob- 
tained from p by tracing out over all spins except spin i 

and where (A) = Tr p- A . The one-tangle r^ ranges 
between and 1. It vanishes for a state such that 
IV') = \4')i '^ I'P) s-nd reaches 1 for maximally entan- 
gled states such as the famous Einstein-Podolsky-Rosen 
(EPR) state H 



iEPR) = i=(in)-uT)). 



(6) 



To measure the two-spin entanglement, we consider 
the concurrence C introduced by Wootters [S^. This 
quantity is computed from the two-spin reduced density 
matrix pij obtained by tracing out p over all spins except 
spins i and j. Next, we introduce the spin- flipped density 
matrix pij — ay®ayp*,jay®ay where p*.j is the complex 
conjugate of p^.j. The concurrence C is then defined by 



Ci_j = max{0, Pl 



M2 



M3 - ^4} = max{0, C* }, (7) 



where the pk 's are the square roots of the four real eigen- 
values of Pi J Pi J and where pi > pi+i- The concurrence 
vanishes, for example, for any state j?/;) = \x)i(^\'P)j'^\'P) 
and reaches its maximum value C — 1 for states such as 
I EPR). 

Let us also mention that r and C are related through 
the Coffman-Kundu-Wootters (CKW) conjecture [43 
stating that 






(8) 



IV. GROUND-STATE ENTANGLEMENT 

We recall in this section the main features of the 
ground state entanglement and we refer the reader to 
Refs. [SJiy^lSal for a detailed discussion. 

As explained in Sec. m the ground state of H is well 
described by a mean-field approach showing that in the 
thermodynamical limit the ground state is given by a 
product state [S^l for which t = C = 0. However, in the 
symmetric phase (A > and ft, < A), we underline that 
the ground state is twofold degenerate in the thermody- 
namical limit whereas at finite-iV, it is unique (except 
for h = 0). In the degenerate case, one must consider 
the thermal density matrix at zero temperature |3J| to 
compute the entanglement properties. This implies in 
particular that the one-tangle in this phase is given by 
T = 1— (/i/A)^, and that C = in the large- A'' limit. How- 
ever, the finite N corrections to this separable asymptotic 
form allow one to capture a nontrivial behavior of the 
rescaled concurrence Cr = {N — 1)C. This rescaling is 
actually the coordination number of each site and takes 
into account the fact that the two-site entanglement, as 
measured by the concurrence, is, in this model, equally 
shared between all sites. 

For a ferromagnetic coupling (A > 0), Cr displays a 
cusplike behavior at the critical point he = X (see Fig. 
^ where Cr goes to 1 in the thermodynamical limit. 

In the antiferromagnetic case, Cr is a smoothly de- 
creasing function of h but is discontinuous at ft. = 
where a first-order transition occurs (see Fig. ^. For 
ft = 0^, Cr reaches 1 in the thermodynamical limit, but 
for ft = 0, the ground state is given by zero total spin 
states for which Cr = 0. 

Except for this very special point, the ground state 
always lies in the maximum spin sector S = N/2 for any 
A. Moreover, at finite N and in the large-field limit, the 
ground state is the fully polarized state in the z direction. 



\%) = \Z)=6, 

for which one also has Cr — 0. 



^W 



T), 



(10) 



Cr 




FIG. 1: Rescaled concurrence of the ground state for A = +1 



(thin line) and A = — 1 (thick hue) {N • 
cases, Cr = for h = 0. 



W^ spins). In both 



In the antiferromagnetic case, the entanglement prop- 
erties of the ground state (restricted to the maximum 
spin sector S — N/2) have also been analyzed by Hines 
et al. 40] through the Von Neumann entropy E{h) for 
bipartite systems. In BEC language, the two subsystems 
correspond to the two components of the condensate, i. e. 
to the states polarized in the ±a; direction. This quan- 
tity has thus nothing to do with the concurrence and its 
behavior is indeed very different: E vanishes iai h — 
(since \'ipf)) = \X)) whereas Cr is maximum (in this spin 
sector) and E is maximum in the large-field limit whereas 
Cr vanishes. 

More recently, another entropy S{h) has been analyzed 
in the LMG model [47| . This entropy is based on a bipar- 
tite decomposition in two subsytenis of size L and N — L. 
Here, S{h) has been found to be maximum at the critical 
point he in agreement with a recent conjecture proposed 
by Hines et al. ^3|- Moreover, a logarithmic scaling has 
been observed and may reveal an underlying ID confor- 
mal field theory ^ for the LMG model. 



V. QUANTUM DYNAMICS 

In this section, we analyze the quantum dynamics of 
two initial states belonging to the maximum spin sub- 
space S = N/2: 

1. the state \X) fully polarized in the x direction 
which corresponds, in BEC language, to a state 
where the atoms are all localized in one of the two 
modes: 

2. the state \Z) fully polarized in the z direction which 
corresponds, in the BEC language, to a state where 
the atoms are equally distributed between both 
modes (catlike state). 

The quantum dynamics of these two states have al- 
ready been widely studied for this model in the BEC 



context iililliilHSIiliilillll. Here, we focus on 
the entanglement dynamics by studying the time evolu- 
tion of t(/i, t) and CR{h, t), which do not depend on the 
spins kept in the trace operation since states belonging 
to the subspace S = N/2 are invariant under the per- 
mutation group. This implies in particular that r and 
Cr, only depend on the magnetization (Sa) and on the 
correlations functions (SaSfs) 55j. 

A. Time scales 

As explained in Ref . |56| , the spectrum of H remains 
discrete in the thermodynamical limit; i. e. ,the mean 
level spacing remains finite when N goes to infinity. For 
equidistant levels, we would expect to have a Poincare 
time that does not depend on N. However, it turns out 
that the Poincare time scales here linearly with N. As 
noticed by Milburn et al. |5J|, the semiclassical dynam- 
ics (see below) mimicks very well the quantum dynam- 
ics up to a given time scale tsc which is of the order 
of the Poincare time. For times smaller than tsc-, the 
system thus remains nonentangled and in the large-TV 
limit, it is described by the classical equations of motion 
(see Appendix IB)| . For times larger than tgc quantum ef- 
fects become important, some revivals of the wavepacket 
are observed, and the latter approach fails. It is thus 
in this regime that interesting and nontrivial entangle- 
ment properties must be analyzed. Consequently, we al- 
ways consider a dimensionless rescaled time fi = A\t/N 
which allows us to investigate the dynamics and to easily 
compare the results for various N . Of course, we have 
checked in each situation that the time-averaged results 
were weakly sensitive to the the time range chosen pro- 
vided this latter was sufficiently large. 

We must also mention recent studies JS^, 113 which 
reveal multiple time scales in the LMG model (called 
the boson-Hubbard dimer model in these papers). How- 
ever, the range of validity of these results is confined to 
h <^ X/N (small-tunncUing-amplitude regime in boson 
language) which is not the regime considered here. Yet 
we cannot exclude larger characteristic time scales but 
such a study is beyond the scope of the present paper 
and we shall leave this question apart in the following. 



B. The \X) state 

It is clear that the entanglement properties of the \X) 
state are insensitive to the transformation h ^^ —h. 
Thus, as explained in Appendix IXI one can, in this case, 
restrict our analysis to A > and h > 0. 

• For h — 0, the \X) state is an eigenstate of H so 
that Tx{h = 0) = and Cx{h = 0) =0 for all time 
t. 

• For X — 0, H is completely separable and thus both 
Tx and Cx vanish. However, for this state, this 



limit is very singular since, for any finite coupling 
A, the large- field limit leads to nontrivial entangle- 
ment properties {tx ^ 1) as can be seen in Fig. 
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FIG. 2: Beliavior of tx as a function of p = AXt/N for differ- 
ent values of the field h = 0.5, 0.55, 0.8, 1 from bottom to top 
{\^+l and AT = 10^). 
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FIG. 3: Same as Fig. Olfor h = 0.2, 0.3, 0.4, 0.45 from bottom 
to top. 

Apart from these limiting cases, we have been unable 
to get analytical expressions of the one-tangle for this 
state. We have displayed in Figs. |21and|31thc behavior of 
Tx for different values of the magnetic field as a function 
of /x. 

For /i <C A, the characteristic energy scale driving the 
dynamics is \/N and hence r is almost periodic in /i 
with period 27r [Fig. |31 (bottom)] . In this regime, quan- 
tum fluctuations are weak and r is a smooth function 
of ^. When ^ increases, the depth of the hollows which 
corresponds to a revival of the wave packet decreases and 
eventually goes to zero in the large /i limit so that, in this 
limit, Tx is almost constant. When h increases, quantum 
fluctuations become more and more important and tx 
becomes more and more "noisy." The maximum of these 
fluctuations is reached for h = A/2 which corresponds to 



the critical field of the semi-classical dynamics |41|, |5^ . 
Below (above) this critical field, {Sx) (which represents 
the difference of population between the two wells in the 
BEC problem) oscillates around a nonvanishing (vanish- 
ing) value. In the other limit h^ X where the magnetic 
field dominates tx goes to 1 and displays a pseudo pe- 
riodic behavior in fi with period At: as can already be 
inferred ioi h — X in Fig . [3 (top). At large times, 
the depth of the cusplike behavior present for /x = 0[47r] 
decreases and r is also almost constant. To go beyond 
this qualitative description, we have computed the time- 
averaged value of the one-tangle tx as a function of h. 
As shown in Fig. 0] the large N limit of tx can be par- 
tially captured by a simple semiclassical approach, es- 
pecially, when quantum fluctuations are weak; i. e., far 
from h = A/2. 
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FIG. 4: Time-averaged value of tx as a function of the mag- 
netic field hiox N = 800, 1600 (A = -1-1). The solid line is 
obtained from the "(semi) classical time-averaged one-tangle" 
as explained in the text. 



For a detailed discussion of this approach, we refer the 
reader to Refs. ^,^J. The main idea is to compute the 
solution of the semiclassical equations for the quantity 
Sa ■ For an initial state polarized in the x direction, these 
equations can be exactly solved following Eilbeck et al. 
and one gets 



sT{t) = 


f cn(2M|fc2), 


(11) 


Syii) = 


-y dn(2M|fc2) sn(2/it|fc2), 


(12) 


sf{t) = 


^^ sn^ {2ht\k^) , 


(13) 



where en, dn, and sn are the Jacobi elliptic functions and 
where fc — X/2h (see AppendixiBlfor details). These S^'s 
are periodic function with a period 2K{k^)/h for fc^ < 1 
and 2K(k^)/X for k^ > 1 where K denotes complete 
elliptic integral of first kind. Note that the case fc^ = 1 is 
singular since the upper and lower limits do not coincide. 



Of course, if we compute the one-tangle using these 
semiclassical expressions, we triviahy find 
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(14) 



^=i-4e(^")^-i-4e(^«^ 



(15) 



P* ■ 

and, consequently, r|? = for all h. However, as said ^-^x.r 

above, the fluctuations of tx are pretty weak at least 

away from h = A/2. Thus, assuming that tx is constant, 

we replace {SaY by (S^) : 



IJJbpJtdt^^ilq^^ 
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Obviously, this crude approximation is to be justified 
only for constant {Sa) and S^, provided the semiclas- 
sical description is meaningful; i. e., in the large- TV 
limit. After time-averaging the S'^'^'s over one period, 
we thus get a "semiclassical time-averaged one-tangle." 
This quantity is depicted in Fig. ^ and as expected it 
is in excellent agreement away from h = X/2. Near this 
point, since the quantum fluctuations are very strong, 
the approximation (|15|l fails and the discrepancy is 
large. In the large-field limit, one has S^^ = for all a 
and thus tx = 1- 

Let us now discuss the concurrence dynamics for the 
\X) state whose behavior is depicted in Figs. [H and 
for different values of the magnetic field. At small field, 
[h < A/2 (see Fig. EJ], Cx r is most often negative 
but displays some peaks (where Cx ^ > 0) which co- 
incides with the hollows in tx and whose amplitude goes 
to zero when /i increases. As a result, at large times, 
the rescaled concurrence Cx,b. always vanishes. When h 
increases, this phenomenon is amplified and, as for the 
one-tangle, Cx n fiuctuates more and more until h — A/2 
where it reaches its minimum value. Above this critical 
field, Cx R increases when h increases, but remains al- 
ways negative at large times. Note that at short times 
there is, for any h, a "peak of concurrence" for which 
Cx n is positive but as time goes on, it is always nega- 
tive. 

To summarize, at sufficiently large times and for any 
values of the magnetic field /i, the rescaled concurrence 
Cx,R always vanishes. 



FIG. 5: Behavior of Cx,r as a function of /^ = 4:\t/N for dif- 
ferent values of the field h = 0.3, 0.4, 0.45 from top to bottom 
(A = 4-1 and A'^ = 10''). The grey line corresponds to the 
threshold C = 0. 
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FIG. 6: Same as Fig. Elfor h = 0.5,0.8, 1,2 from bottom to 
top. 



in the dynamics of \Z). Thus, to investigate all possible 
situations, we consider the case A > for all h, the region 
A < being obtained by transformations discussed in 
Appendix 1X1 

• For h — 0, H is trivially solvable as well as the time 
evolution of the \Z) state. As detailed in Ref. [55|. 
one finds, in this case 



C. The \Z) state 

The \Z) state is an eigenstate of the operator Yii'^l 
which commutes with H. As a consequence, at all times, 
one has (Sx) = {Sy) = 0, so that 



^^ = l-]^('5. 



(16) 



In addition and contrary to the \X) state, it is obvious 
that the direction of the magnetic field h plays a role 



Tz{h^O,t) 
CzMih^O,t) 



1-cos(m/2)2(^-i), 



- cos(^)---l 



([cos(Ai) 



^-2 -if 



(17) 



(18) 



[4cos(Ai/2)^-2gijj(^/2)]^^^ ^ 



where, as previously, fi — 4Xt/N. The expression 
for Cz,R{h = 0, t) given above is obtained by noting 



that the condition 



24^) 



< 



pIMI^ 
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I (2)1 

\Pl4 



(19) 



is always satisfied. Thus in the large- iV limit, tz 
is equal to 1 except for fi — [27r] where tz — 
and, as expected, the concurrence vanishes for all 
fi. However, if one considers the rescaled concur- 
rence, a close inspection of expression H19|l shows 
that a periodic pulse also occurs for /i = [27r] 
whose width becomes sharper and shaper as N in- 
creases and whose maximum height goes to 1 in the 
large- TV limit. 

At this stage, it is interesting to note that tz and 
Cz.R are strongly correlated in the sense that, 
when the one-spin entanglement is large (and here 
maximum), the two-spin entanglement is small 
except for /i ~ [2tt] . This is a relatively surprising 
result which is also encountered in other regimes 
(see below). Indeed, usually, the trace operation 
destroys some correlations and one expects that 
the absence of two-spin entanglement might not be 
able to generate one-spin entanglement. Neverthe- 
less, as discussed in Sec. IIIII this phenomenon is 
met for the |GHZ) state, but here it is encountered 
for a nonstationary state which preserves this 
feature as time goes on. 



For A = 0, the \Z) state is an eigenstate of H so 
that rz(A = 0) = and Cz(A = 0) = for all time 
t. Contrarily to the \X) state, the limit A = coin- 
cides with the high field regime at finite A. Indeed, 
as discussed in Ref. Q , the | Z) state is the ground 
state of H for h > X in the thermodynamical limit, 
so that Tz{h > A) ~ and Cz{h > A) ~ 0. 

However, the rescaled concurrence has a nontrivial 
behavior that can be captured due to its relation 
with the spin squeezing parameter jSSj 



e 



N 



(20) 



which measures the spin fluctuations of a correlated 
quantum state |58| . The subscript n± refers to an 
axis perpendicular to the mean spin (S) where the 
minimal value of the variance is obtained. 

Indeed, as shown in Refs. [5^ iGfl, for any state 
belonging to the subspace S — N/2 and eigenstate 
of the operator JJ- (j\ , one has 



e 



l-C 



R, 



(21) 



provided the reduced density matrix elements of 
the state considered verify |pi4 | > P22 ■ We have 
checked that actually this inequality is always sat- 
isfied for any h and for all the times under investi- 
gation. 



• The relation H21|l thus allows us to deduce the 
rescaled concurrence using the frozen-spin approx- 



imation used in Ref. 
different regimes. 



|49j | to compute ^ in two 



For /i > A, as explained above, the \Z) state is 
the ground state of H in the infinite N limit. 
So, even at finite iV, one expects (Sz) to be 
almost constant in this region. The so-called 
frozen-spin approximation precisely consists 
in replacing, in the Heisenberg equation, SJt) 
by its initial value N/2. Then, following [i^ 
and using Eq. H21I) , one gets an analytical ex- 
pression for the rescaled concurrence, 



C. 



Z,R 



1 



o , . (JO • 2 / \ 

cos (ijjt) + ——- sm [Lot) 



(22) 



with cj = 2y^h{h — A). 



— For /i < 0, the situation iSjin fact, similar. 
Indeed, as shown in Ref. j^, the \Z) state 
is the ground state of H{X < 0, /i > 0) in the 
thermodynamical limit. As detailed in Ap- 
pendix^ this implies that it is the highest- 
energy state of H{X > 0,h < 0). Thus, the 
frozen-spin approximation can also be used in 
the region {h < 0) and, according to [43, one 
then obtains 



C 



Z,R 



1 



COS [Lut) H sm (ujt) 



(23) 



Of course, this approximation is more and more 
valid in these regions {h < and h > X) when 
N increases. Conversely, at fixed N, it is safe for 
either /i > A or /i < 0. 

The only region where the | Z) state is not an eigenstate 
of H in the thermodynamical limit is < ft- < A. We 
have plotted in Fig. [7|the behavior of tz as a function of 
the rescaled dimensionless time ^ = AXt/N for different 
values of the magnetic field h in this range. 

As explained above, r^ ~ for /i < and h > X since, 
there, \Z) is an eigenstate of H in the large-A^ limit. By 
contrast, for h G [0, A], tz is an oscillating function of /i 
with an average value that is nonvanishing as can be seen 
in Fig. [3 The time-averaged value of tz for different N 
is shown in Fig. |S1 

As can be seen, rz is a decreasing function of N at fixed 
h except for /i = as discussed above. To determine its 
large- A'^ behavior, we have computed tz for several h, as 
a function of N, up to A^ = 5000. The results displayed 
in Fig. I^do not allow us to extract an asymptotic value 
but it is likely that for h £ [0, A[, limjv_oo ^^ 7^ 0. 

As for the |X) state, we could have tried to reproduce 
this behavior by a classical analysis, but unfortunately, 
the \Z) state is a fixed point of the classical equations 
of motion (IB8II and IBQI). Thus, one must in this case 
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FIG. 7: Behavior of tz as a function of /i = AXt/N for /i 
0.01 (grey) and h = 0.9 (black) (A = +1 and A^ = 10^). 
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FIG. 9: Behavior of tz as a function of A'" for h = (top) to 
1 (bottom) with a step Aft = 10"^ (A = +1). 
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FIG. 8: Time-averaged value of tz as a function of the mag 
netic field ft for iV = 400, 800, 1600 (A = +1). 



use a real semiclassical treatment analogous to that pre- 
sented in |42| which consists in averaging over classical 
trajectories near this fixed point. However, as shown in 
Fig. [7| T^ is a strongly oscillating function so that the 
approximation (|15|l is not valid and we did not compute 
the "semiclassical time-averaged one-tangle" for the \Z) 
state. 

We now focus on the concurrence dynamics for the 
\Z) state, whose behavior is shown in Fig. ^|for h = 
5 A. As can be seen, the frozen-spin approximation is, in 
this case, in good agreement with the exact (numerical) 
results. For h G [0,A], C^^ is always negative so that 
Cz.R = 0. It is worth noting that there is once again an 
"anticorrelation" between the one-tangle and the rescaled 
concurrence since tx ^ when Cz,b. = for /i e [0, A], 
whereas for h g] — oo,0] U [A, oo[, one has tx — and 
Cz,R 7^ 0. Indeed, in this latter region where \Z) becomes 
an eigenstate in the thermodynamical limit, one always 
has a non vanishing rescaled concurrence as can be seen 
in Fig. 111! where we have displayed the time-averaged 
value of Cz,R as a function of h. At large \h\ the frozen- 
spin approximation allows us to extract the asymptotic 



'Z,R 




FIG. 10: Behavior of Cz,r as a function of /j, = AXt/N for h = 
5 (A = -1-1 and A'^ — 10^). Dots are obtained from the exact 
numerical diagonalizations and the solid line is computed with 
the frozen-spin approximation 12211 . 



behavior Cz,b. ~ |A//i|. 



VI. CONCLUSION AND OUTLOOKS 

We have analyzed the one-spin and two-spin entan- 
glement dynamics for two different states. In the case 
of a state fully polarized in the x direction, we have 
found a nontrivial behavior of the one-tangle whose time- 
averaged value increases monotonically when the mag- 
netic field increases whereas it has a vanishing concur- 
rence for all h. By contrast, for an initial state fully po- 
larized in the z direction, the one-tangle is nonvanishing 
only for h Cz [0, A] where its concurrence vanishes. Apart 
from this region, the one-tangle vanishes in the thermo- 
dynamical limit (since \Z) is then an eigenstate of H) but 
the rescaled concurrence has a nontrivial time-averaged 
value which can be understood using the frozen-spin ap- 
proximation introduced in Ref. |42j . 

Contrary to the ground-state entanglement, we did not 
find a nontrivial singularity at the critical point h = X 




FIG. 11: Time-averaged value of Cz.r as a function of the 
magnetic field /i for Af = 800, 1600 (A = +1). The solid 
line corresponds to both branches (/i < and ft > 1) of the 
frozen-spin approximation (FSA). 



but only some enhanced fluctuations near h = A/2 for 
the \X) state. For the \Z) state, we have not detected 
any peculiarity at this point except that dhJz seems to 
be minimum there whereas dhJx is rather maximum. 

We wish to emphasize that some experiments may be 
possible to investigate these dynamical entanglement, es- 
pecially the concurrence. Indeed, as explained in the 
Sec. IV CI the rescaled concurrence can be related to the 
spin squeezing parameter as established in Refs. 59, 60j. 
Thus, the behavior of the Cr as a function of the field 
should, in principle, be measured for both the ground 
state and \Z{t)). Nevertheless, such an exciting perspec- 
tive is of course conditioned by the possibility to measure 
the squeezing parameter. Bose-Einstein condensates are 
certainly the best candidates for such experiments since 
it has already been possible to measure the angular mo- 
mentum in such systems |61|. 



Ackno'wledginents 

We would like thank B. Dougot, J. Dukelsky, R. 
Mosseri, and D. Mouhanna for fruitful discussions. 



it is then straightforward to show that 

2"-l 
i=0 



satisfies 



H{\,-h)\^) = E\^). 



(A3) 



(A4) 
(A5) 



Here, we have used a natural coding of the Hilbcrt space 
states where the number i (and its associated state |z)) 
is in a one-to-one correspondence with its binary decom- 
position expressed in terms of "f (1) or [ (0) spins. These 
identities show that the spectrum of H is odd under the 
transformation A —t —A and even under h — > —h. More- 
over, they also imply that 



H{-\^h)\<t>) 



(A6) 



so that i/(+A, +/i) and H{—\, —h) have the same eigen- 
states but with opposite eigenenergies. 

Let us now consider a state {ip) whose spectral decom- 
position on the eigenstates of H reads 



2"-l 
3=0 



(A7) 



where all the a^'s are real numbers, and an observable A 
satisfying 

(</.,v|^|0,.) = (0,1^1^,,), (A8) 

for all {j,j')- Then from Eq. HA6() one has 

ii^+MtUl^+At)) - {i^-At)\AH-At)). (A9) 

where \i{j±,±{t)) = e"''-^(^^'=^'')*|i/'). In addition, it is 
clear that the density matrix matrices p+,+ and p-,- are 
identical. 

APPENDIX B: SEMICLASSICAL APPROACH 



APPENDIX A: SPECTRUM SYMMETRIES AND 
THEIR IMPLICATIONS FOR THE DYNAMICS 

Let us consider an eigenstate of the Hamiltonian H, 

2~-l 



i=0 



such that 



HiX,hM^E\^) 



(A2) 



The Heisenberg equations for the spin operators are 
given by 



ijx — z^nijy 



2A, 



by — Ztlbx H "Tfybzbx + bxb. 



N 



2A 



i^z — 'ZTK^y^x ~\~ Jxby)- 



(Bl) 

(B2) 

(B3) 



These equations cannot be solved exactly for any (A, h) 
and any initial conditions. To investigate the spin dy- 
namics, we follow Milburn et al. ISJl and consider the 



semiclassical dynamics (mean-field approximation) asso- 
ciated with these equations by considering the "double- 
well" parametrization (Schwinger-like) for the spin vari- 
ables which reads 



f(N 



\bi\ 



iN 



(6J62 - hb\ 



'2) J 



N 



{b\b2+b^bl). 



(B4) 
(B5) 
(B6) 



The complex variables &i and 62 satisfy the constraint 
|6i|2 + |62|2 = l and obey 



bj = ihb^-j + 2iX\bj\'^bj, 



(B7) 



for j = 1,2. This type of equation has been studied in 
detail by Eilbeck et al. so that we only sketch here the 
main lines of the solutions and refer the reader to Ref. 
[531 for details. Setting bj = Oj-e'^^ (for j = 1,2) and 
6 = 9i — 02, the latter nonlinear equation system HB7|I is 
rewritten as 



2\l/2 




Then, using the conservation of energy, 

E{ai,e) = -X[l-2al{l-al)] 



-2hai (1 



,2\l/2 



008 6*, 



1 

N 



Hiai 



-v.-^ 



(B8) 
(B9) 

(BIO) 
(Bll) 



we can eliminate 6 to get a closed equation for ai. 

From now on, we consider the initial conditions ai (0) = 
1 and 6{0) = which describe a state polarized in the 
—X direction; i.e., the classical counterpart of the | — X) 
state. This state is not the one considered in the Sec. 
IV Bl since it lies in the opposite direction, but for clarity, 
we have used the same notation as in Refs. |5J,|53- One 
thus has 

E{ai,e) = -X, (B12) 

so that from the energy conservation ljB10|l . one gets 

cos6'= -ai(l-a?)i/2. (B13) 

Next, using Eq. I)B8|I and setting al — (1 + y)/2, it is 
straightforward to show that 



^1 =1 
2h 



k^ + (2fc2 - l)y2 „ fc2^4^ 



(B14) 



where k — X/{2h). The solution of this equation is the 
Jacobi elliptic function 



y{t) = cn(2M|fc2 



(B15) 



which yields 

S^t) = y(l-2a2)=_^cn(2/ii|P), (B16) 

Sy{t) = -NaiCLi/h = — dn{2ht\k^) sn{2ht\k^) , (B17) 

S^_{t) = ^alil - aif ^ ^ sn^2ht\k^) . (B18) 

Finally, since we are not interested in the | — X) but in 
the \X) state, we just have to change the sign correctly 
and we get the expression given in Sec. IV Bl 



APPENDIX C: ONE-SPIN AND TWO-SPIN 

REDUCED DENSITY MATRICES OF FULLY 

SYMMETRIC STATES 

Let us consider a state 



+N/2 

E 

M=-N/2 



aM\N/2,M), 



(CI) 



where {\S, AI)} is an eigenbasis of S^ and Sz, and its den- 
sity matrix p = lip) {ip\. Here, we restrict our discussion 
to states belonging to the maximum spin sector S = N/2 
which are relevant for BEC. 

In this subspace, all states are invariant under the 
permutation group Sn, so that one can easily compute 

the matrix elements of the reduced density matrices p] 

(2) 
and pI ' . Note that the permutation symmetry implies 

that pI does not depend on i, and that pl ■ does not 
depend on i and j, so that wc omit these indices in the 
following. 

As is well known, the one-spin reduced density matrix 
is easily expressed in terms of the one-spin correlation 
functions or, more precisely, in terms of the (o'q.)'s. In the 
eigenbasis {| t), | i)} oiaz, a straightforward calculations 
gives 



p['' 



,M 



M 

(1) 1 Y^ I 12*-^ 

P22 = ^-I^I^mI ^, 

M 

(1) V- . J{N-2M + 2){N + 2M) 

P12 = Z^ctMa-M+i 



M 



2N 



(C2) 
(C3) 

,(C4) 



•+u (1) (1)* 

withp)^ ^py ■ 

The two-spin reduced density matrix can also be easily 
expressed in terms of the two-spin correlation functions 
(see Ref. 113). In the eigenbasis {| tt), I Ti), I it), I ii)}, 
one has 



(2 



Pi 



i^ = E 



„ (iV + 2Af)(JV + 2M-2) 



M 



4:N{N-1) 



(2) 
Pl2 



(2) 
Pl4 



(2) 
P24 



^«Afa;,+i(iV + 2M-2) 



M 



y^{N-2M + 2){N + 2M) 

^ miN-i) ' 

Y, aMaM+2 V{N + 2M)iN + 2M~2 



(C6) 



M 



y/{N - 2M + 2){N - 2M + 4) 
J2aMa*M+iiN-2M) 

M 



(C7) 



(2) 
P22 



(2) 
P44 



V^(A^ - 2M + 2){N + 2M) 



E 

M 

E 

M 



aM 



aM 



4iV(iV-l) 

^{N ~ 2M){N + 2M) 
AN{N - 1) ' 

,(iV-2M)(iV-2M-2) 



4^(A^-1) 
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(C8) 
(C9) 

(CIO) 



with pf^ = pf2\ P23 

(2) (2)* 



P§ = P221 P34 = P24 > and 
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